Evidence for gapped spin-wave excitations in the frustrated Gd2Sn207 pyrochlore 
antiferromagnet from low-temperature specific heat measurements 
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We have measured the low-temperature specific heat of the geometrically frustrated pyrochlore 
Heisenberg antiferromagnet Gd2Sn207 in zero magnetic field. The specific heat is found to drop 
exponentially below approximately 350 mK. This provides evidence for a gapped spin-wave spectrum 
due to an anisotropy resulting from single ion effects and long-range dipolar interactions. The data 
are well fitted by linear spin-wave theory, ruling out unconventional low-energy magnetic excitations 
in this system, and allowing a determination of the pertinent exchange interactions in this material. 



In the insulating R2M2O7 magnetic pyrochlores, 
the rare-earth ions (R"^^) sit on a lattice of corner- 
sharing tetrahedra, resulting in geometrical frustration 
for antiferromagnetic nearest-neighbor exchange between 
Heisenberg spins. As such, these materials display a 
wide variety of interesting and unusual behaviors. Ex- 
amples of observed phenomenology include spin ice [if, 
spin glass 0, and spin liquid [ll behaviors as well 
as long-range magnetic order 0,13, H, [ol. [Tol . \v\\ and 
perhaps some novel type of hidden order [12] ■ A com- 
mon trend observed in all of these systems is an ap- 
parent persistence of spin dynamics down to the low- 
est ternperatures as revealed by muon spin relaxation 
(AtSR) [a,S,l3,ll,[i,[illJll, MSssbauer spectroscopy ^ 
and neutron spin echo 1] experiments. The origin of per- 
sistent spin dynamics (PSDs) is a major open question in 
the study of highly frustrated insulating magnetic oxide 
materials Q. 

Perhaps what is most perplexing, is that PSDs have 
been found in pyrochlores, such as Gd2Ti207 (GTO) [9, 
m and Gd2Sn207 (GSO) 0, HJiIL which, according to 
neutron scattering experiments [6!,[l0'|, display long-range 
magnetic order below a critical temperature Tc — 0.74 K 
for GTO [6] and = 1.0 K for GSO [13]. This is 
highly unusual since conventional wisdom suggests that 
collective magnon-like excitations, hence spin dynamics, 
should freeze out in the limit of zero temperature. Sim- 
ilarly, the Gd3Ga50i2 antiferromagnetic garnet (GGG), 
which displays rather extended magnetic correlations be- 
low 140 mK [3] that are now theoretically rational- 
ized pj], also exhibits PSDs below 100 mK [11]. It is 
tempting to speculate that the unusual low-energy ex- 
citations giving rise to PSDs may result from a rem- 
nant of the high frustration that these systems possess 
and the lack of long-range order that would occur were 
it not for perturbativel y sm all exchange interactions be- 
yond nearest-neighbor [l3] , dipolar interactions [3, [13] 



or single- ion anisotropy [20]. 

We see Gd2Sn207 (GSO) as a key system to inves- 
tigate in seeking to convincingly demonstrate, via mea- 
surements of the low temperature magnetic specific heat, 
Cm{T), the existence of unconventional low-energy exci- 
tations, as suggested by the observation of PSDs. Previ- 
ous Cm measurements on this material between 350 mK 
and 800 mK (see Fig. 1) were found to be parametrized 
by a Cm{T) - law Such a temperature de- 

pendence of Cm is unusual since conventional antifer- 
romagnetic magnon excitations without a gap lead to 
Cm{T) ~ or, with an anisotropy energy gap A, to 
Cm{T) ~ exp(— A/r). Such an unconventional Cm{T) 
behavior further argues for the existence of unusual low- 
energy excitations in GSO. Interestingly, the related 
GTO, which also exhibits PSDs, displays Cm{T) ~ 
for 100 mK < T < 500 mK [§]. 

In contrast, a recent theoretical work [l^ argues 
that the lowest temperature ~ 350 mK considered in 
Ref. [H] for GSO corresponds to the upper tempera- 
ture limit above which multi-magnon excitations start 
to proliferate, obscuring the genuine single-magnon/low- 
temperature regime. In particular, calculations using the 
microscopic spin Hamiltonian 1^, 18, 1^ that describes 
the experimentally observed ground state of GSO [13] 
predict that Cm should begin to drop away exponentially 
exactly at or just below 350 mK as gapped magnetic ex- 
citations become quenched at lower temperature |19i] . 

In this letter we present an investigation of the na- 
ture of the low-energy spin excitations in GSO through 
measurements of Cm(T) down to 115 mK. The results 
show an exponentially dropping specific heat at low tem- 
perature and, when compared with theoretical calcula- 
tions [13], confirm the picture of conventional gapped 
antiferromagnetic magnons in Gd2Sn207. 

PolycrystaUine samples of Gd2Sn207 were prepared by 
solid state reaction as in Ref. 2l|. Gd203 and Sn02 
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were mixed in stoichiometric ratio, pressed into pellets, 
and then fired in air for 48 hours at 1350°C. The pellets 
were then reground and the process repeated. Powder 
X-ray diffraction spectra taken after the two firings were 
indistinguishable. All peaks were indexed to space group 
Fd3m, with the size a of the cubic unit cell a = 10.460 A, 
with no impurity peaks detectable at the 1% level. High 
temperature susceptibility measurements yield a mag- 
netic moment of 7.95/iB for Gd'^^ in GSO, close to the 
expected free-ion moment of 7.94/is. 

Specific heat measurements were performed using the 
quasi-adiabatic technique. A 1 kfl Ru02 resistor was 
used as thermometer and calibrated to a commercially 
calibrated germanium resistance thermometer. A 10 kfi 
metal-film resistor was used as a heater. The sample was 
suspended from thin nylon threads and the heater and 
thermometer were fixed directly to the sample. Leads to 
the thermometer and heater were made from 6 ^m di- 
ameter, 1 cm long NbTi wires. The sample was weakly 
heat-sunk to the mixing chamber of a ■^He/'*He dilution 
refrigerator with Pto.92Wo.08 wire, chosen for its insignifi- 
cant contribution to the addendum. The addendum, due 
to other components such as heater and thermometer, 
was determined to be less than 2% of the sample's heat 
capacity in the worst case at around 120 mK and thus 
does not effect the results of our analysis. 

Slow thermal relaxation within the sample was ob- 
served (with a relaxation time constant ri ^ 120 s at 200 
mK). The thermal link was chosen so that the relaxation 
time of the sample temperature to cryostat temperature, 
T2, would be much longer than ri, thereby minimizing 
the temperature gradient within the sample. A double 
exponential form was fit to the sample's temperature as 
a function of time after a heat pulse and the longer expo- 
nential was extrapolated to zero time, giving a measure 
of the sample's heat capacity. The data presented here 
was taken using a long time constant T2 ~ 20ri . Another 
experiment performed with T2 ~ 5ti resulted in a slightly 
noisier measurement. Yet, the specific heat in both mea- 
surements overlapped within the estimated error bars of 
the data, ruling out significant inaccuracies due to the 
slow thermal relaxation in the sample. 

Our specific heat results, shown in Fig. [1] are largely 
consistent with previous work though there is a 10 to 
15% systematic discrepancy between the data sets. Our 
measurements, however, extend to considerably lower 
temperature 115 mK) allowing us to test the proposal 
of gapped magnon excitations [l9|. We observe below 
~ 350 mK a deviation from the behavior describing 
the data between 350 mK and 800 mK as previously re- 
ported 0. The specific heat decreases also faster than 
the power-law expected for conventional gapless anti- 
ferromagnetic magnons. This already suggests that the 
specific heat may be dropping out exponentially, indicat- 
ing a gapped spin-wave spectrum. 

At the lowest temperatures an upturn in C(T) be- 
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FIG. 1: (Color Online) Specific heat of Gd2Sn207 as a func- 
tion of temperature from this work (blue circles) between 115 
and 800 mK. Previously measured data at higher tempera- 
tures(green diamonds) the power law previously pro- 
posed (dotted red line) and a power law (dashed blue line) 
are also plotted. The upturn seen below 150 mK results from 
the nuclear electric quadrupole interaction (solid line). 



comes apparent and can be ascribed to a nuclear contri- 
bution, C]\[(T), stemming largely from the nuclear elec- 
tric quadrupole interaction and, to a much lesser ex- 
tent, to the nuclear hyperfine interaction. By properly 
parametrizing this term and subtracting it from the to- 
tal specific heat, we can isolate the contribution from the 
electronic moments. ^^^Gd Mossbauer spectroscopy ex- 
periments find a nuclear quadrupole electric splitting of 
-4.0 mm s~^ 13, 2^- More specifically, the four nuclear 
states are split with energies 0, 0.05, 12.1 and 15.9 mK as 
a result of both nuclear quadrupole and nuclear hyperfine 
interactions [l3| . The only other naturally occurring iso- 
tope that contributes to the nuclear specific heat is ^^''Gd 
which, like ^^^Gd, has spin I = 3/2. The ratio of the nu- 
clear dipole moments is fij^'^ / n]^^ ~ 0.8 and the ratio 
of quadrupole moments is Q]^'^/Q}^^ ^ 1-1 ^23]. The 



quadrupole interaction is expected to be proportional to 
//? where the local [111] ^-direction points to the cen- 



and ly terms 



ter of the Gd tetrahedra. There are no /; 
in the nuclear electric quadrupole interaction due to ax- 
ial symmetry of the Gd magnetic site [S^l • The nuclear 
hyperfine interaction is proportional to • J- How- 
ever, in the Palmer-Chalker ground state of GSO [isj . 
the electronic moments are ordered perpendicular to the 
local z-direction, so we can simply write the nuclear hy- 
perfine term proportional to (InIx- Thus the total nu- 
clear Hamiltonian is Ti.N — cl^^nI-x + bQ^I^. We have 
calculated the nuclear contribution, Cn{T), to C{T) by 
choosing coefficients a and b such that the resulting eigen- 
values of T^AT match the above energy splittings for ^^^Gd. 
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FIG. 

The low temperature limit of the magnetic specific 
heat is parametrized by a phenomenological form 
Cm oc (l/r2)e-^/'^, which appears as a straight line. 
The specific heat is plotted before (squares) and after 
(circles) subtraction of a nuclear contribution. The linear fit 
is performed between the lowest temperature measured and 
250 mK and gives A = 1.40 ± 0.01 K. 
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FIG. 3: (Color online) Fit predicted by linear spin-wave the- 
ory to Cm{T) with resulting values of the exchange interac- 
tions. The best fit gives J2/IJ1I = 0.01 and J31/IJ1I = -0.01. 
Deviation from theory occurs close to the transition, Tc — 
1.0 K, where magnon-magnon interactions become important. 



The interactions are then scaled by the ratios of nuclear 
moments to get the contribution from ^^^Gd. The calcu- 
lated specific heats are multiplied by the relative natural 
isotopic abundances and added to give the total Cn{T). 

At temperatures much higher than the nuclear energy 
splittings, the nuclear contribution behaves as Cjv = 
AT~^ with, based on the calculation above, A ~ \.?>^ x 
10"'' JK/mol Gd. However, because of uncertainties 
in the nuclear moments, we keep ^ as a free parame- 
ter. Cn{T) can be subtracted from the total C{T) at 
these temperatures since the hyperfine coupling is only 
a few mK. As a first rough indicator of a gapped exci- 
tation spectrum, we note that the electronic part of C 
can be parametrized at low temperature by Cm(T) oc 
{1/T^)e~'^^'^ as in conventional colinear antiferromag- 
nets with single ion anisotropy. A plot of log(CT^) ver- 
sus 1/T should then give a straight line and this appears 
to be the case if A is chosen to be 1.63 x 10^'* JK/mol Gd 
as shown in Fig. [2l which is close to the value calculated 
above. The discrepancy is reasonable considering the un- 
certainty on the parameters in the calculation and the 
limited amount of data in the temperature range where 
the nuclear contribution becomes dominant. The phe- 
nomenological "gap" A so obtained is A = 1.40±0.01 K. 
We now go on to describe a microscopic calculation which 
further supports the conclusion that a gapped magnon 
excitation spectrum exists in GSO. 

Gd2Sn207, and the related material Gd2Ti207, are 
well approximated as isotropic, Heisenberg antiferromag- 
nets as the Gd^+ ions have half- filled shells {L = 0, 
S = 7/2). Neutron scattering experiments have identified 
a fc = (0, 0, 0) ordered ground state [l3| - the so-called 



Palmer- Chalker state [13|. Thus Gd2Sn207 is particu- 
larly well suited to the standard protocol of, after hav- 
ing identified the ground state, computing the second- 
quantized (ma gno n) low-lying spin excitations around 
that state [l9|, |2J]. The nearest neighbor exchange in- 
teraction is estimated from the Curie- Weiss constant 
0CW = -8.6 K by Ji = 30cw/zS{S+l) = -0.273 K Q 
where z = 6 is the number of nearest neighbors. The 
strength of the dipolar interaction is derived from the 
size of the magnetic moments and the geometry of the 
lattice 1^, 1^, 24 1 . Though the moments are fairly 
isotropic, there exists a small crystal field anisotropy re- 
sulting from second order admixing from the electronic 
L=0 state to excited manifolds [25|). The crystal field 
parameter i?^ = 47 mK was used in this calculation as 
discussed in Ref. [lit . 

A general two-body spin interaction Hamiltonian 
which includes isotropic magnetic exchange interactions 
up to third nearest neighbor and anisotropy in the form 
of interactions with the local crystal field as well as 
long range dipole-dipole interactions [3, [H, H, 24 1 was 
considered. The techniques of linear spin-wave theory 
Bin were employed to diagonalize the Hamiltonian 
and arrive at a Bose gas of non-interacting magnons cre- 
ated (annihilated) by operators aj^(fc) (a^(fc)). 



4(fc)aa(fc) 
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where Ti^*^^ is the classical ground state energy and ea{k) 
are the spin- wave dispersion energies. The pyrochlore 
lattice can be viewed as a FCC lattice with a 4-site ba- 
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sis (e.g. a Gd tetrahedron). The summation runs over 
all wavevectors k in the first Brillouin zone of the FCC 
lattice and the subscript a runs over the four sublattice 
sites. CmjT) is obtained from standard thermodynamic 
relations US 



C„(T) = fcB/32^^[£a(fc) 



exp[/3ea(fc) 



[exp[/3e,(fc)] - 1] 



(2) 



where /3 = 1/k^T, fee the Boltzmann constant and T the 
temperature. 

Using a maximum likelihood estimator as in Ref. [l^ . 
we have confirmed that the specific heat is fit quan- 
titatively well by a model that includes weak ferro- 
magnetic second nearest neighbor and antiferromagnetic 
third nearest neighbor interactions {J2/\Ji\ = 0.01, 
J3i/\Ji\ = —0.01) as seen in Fig. O Note, however, 
that qualitatively good fits are observed for a substan- 
tial region of the J2 — J31 plane, and there can be lit- 
tle doubt about the presence of a gap to spin excita- 
tions at low temperatures. With the above values of 
J2 and J31, the calculated dispersion relation shows a 
true spin wave gap, Asw ~ 1-24 K, at the F point 
(fc — 0). The difference between this Asw value and 
the A = 1.40 K value above arises from the fact that A 
is obtained on the basis of a phenomenological paramet- 
ric form, Crn{T) oc {1/T^)e~^/'^ , while Agw is a truly 
microscopic gap. With the values of Ji, J2, J31 above, 
the reduction of the classical order parameter (5') = 7/2 
due to quantum fluctuations is (SS) w 0.11 [lol. [23|. 

In conclusion, the experimental results presented here, 
in conjunction with those from a microscopic theoretical 
calculation that builds on a model that characterizes the 
long range ordered ground state of Gd2Sn2 07 fid, H, H, 
241 ■ leave very little doubt that the bulk excitations in 
this system are conventional gapped collective magnons. 
The following question therefore arises: what is the mi- 
croscopic origin of the temperature-independent fiSH re- 
laxation rate Q and the higher than expected hyperfine 
temperature observed with Mossbauer spectroscopy [l3|? 
Such gapped spin- wave excitations would typically be ex- 
pected to lead to a sharply dropping spin relaxation rate. 

The peculiar behavior of the specific heat of GTO 
also remains a mystery for, in contrast to GSO, it does 
not exhibit a gap down to 100 mK as might be 
naively expected. It would be valuable to investigate 
the possibility that the power law in GTO also gives 
way to an exponentially decaying specific heat at lower 
temperature. GSO and GTO are perhaps at this time 
the two pyrochlore materials most amenable to match- 
ing theory with experiment 1^, 3]. Hence determining 



vide insight into the cause of persistent spin dynam- 
ics in other more exotic, highly frustrated systems such 
as Tb2Ti207 0, Tb2Sn207 Yb2Ti207 [H and 

Gd3Ga50i2 (GGG) 
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